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,J^ I Abstract: Recently, several authors studied small quantum systems weakly coupled to free boson or fermion fields 

M-i at positive temperature. All the rigorous approaches we are aware of employ complex deformations of Liouvillians 

(— I or Mourre theory (the infinitesimal version of the former). We present an approach based on polymer expansions 

of statistical mechanics. Despite the fact that our approach is elementary, our results are slightly sharper than those 

contained in the literature up to now. We show that, whenever the small quantum system is known to admit a 

Markov approximation (Pauli master equation aka Lindblad equation) in the weak coupling limit, and the Markov 

approximation is exponentially mixing, then the weakly coupled system approaches a unique invariant state that 

CO is perturbatively close to its Markov approximation. 

> 

o 

o 1 Introduction 

\n 1.1 Motivation 

o ; 

^D . Quantum systems consisting of a small subsystem (say, an atom) and a large component (say, a heat bath) have 

\ I ' received a lot of attention lately, sparked by the elegant results of I231I2T1 l26l . The challenge in this problem 

I> is to prove that the subsystem thermalizes under influence of the heat bath, this property will be called 'return 

to equilibrium', or simply 'RTE' hereafter. The quoted works show that for typical such systems (in a precisely 

defined sense of typicality) the subsystem is close to equilibrium for most times. However, if one wants to study 

C^ in more detail the subsystem dynamics one needs to resort to concrete models. Indeed, such subsystem-reservoir 

models have been successfully and rigorously studied since the late 90's (we refer to fT7ll2| for early results in this 

field) under the assumptions 

1) that the heat bath consists of a free (and hence explicitly solvable) field. 

2) that the coupling between subsystem and heat bath is small compared to the energy scales of the subsystem 
Nevertheless, these results either impose rather strong assumptions on the form of the system-reservoir coupling, 
or they become quite involved technically. The aim of the present paper is to develop an intuitive and simple 
approach for RTE in this case. To explain our result, let us first recall that these systems were already studied in 
the 70's from the point of view of quantum master equations; B. Davies [6] pioneered the rigorous derivation of 
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master equations in this framework, thus making precise earlier heuristic ideas of I. Prigogine and P. Van Hove. 
The master equation is derived, under mild conditions on the form of the coupling, by scaling time t ast ^ X^^ (A 
is the coupling strength) and taking A ^ 0. It exhibits all irreversible phenomena expected in such model systems 
and as such, it has inspired many researchers in open quantum systems. However, it does of course not yield 
information on the long time (longer than A^^) behaviour of the system. 

We prove that, if the condition necessary for the derivation of the master equation is satisfied and the master 
equation is exponentially ergodic (exhibits exponentially fast return to equilibrium), then the system thermalizes 
in the long time limit, for small but nonzero coupling strength A (more generally, it reaches a steady state, since 
we do not assume that the heat bath(s) is(are) in equilibrium). Moreover, we give an explicit bound on the speed 
of convergence towards the steady state. The necessary condition is the integrability in time of certain correlation 
functions of the free field. To our best knowledge, this condition is weaker than that of other RTE- results in the 
literature. 

Acknowledgements. This work was done while W.D.R was a post-doc at the university of Helsinki. We thank the 
European Research Council and the Academy of Finland for financial support. W.D.R. is grateful to Kevin Schnelli 
for serious proofreading and for pointing out several errors in the manuscript. 

1.2 Setup 

Let .^ be a finite-dimensional Hilbert space (modeling the small system) with a Hamiltonian Hs (a Hermitian 
matrix). To describe the field that plays the role of reservoir, we first pick a finite, discrete hypercube A = A^ = 
Z'^n {~L, L]"^ for some L £ N and we enclose the field in the volume A (we could as well choose A to be a box in M'^). 
Since we will mainly use the Fourier transform, we define the set of (quasi-)momenta A* = (7rZ/L)'^n (— tt, tt]'^. The 
dynamics of one reservoir excitation is given by the one-particle dispersion relation oj-^ (q) and the Hamiltonian of 
the whole field in A is given by 

H^ := J2 ^^ilK^, (1-1) 

acting on the bosonic (symmetric) Fock space J^ = r(^^(A*)). Here a*aq are the creation/annihilation operators 
of a mode with (quasi-)momentum q e A*, satisfying the canonical commutation relations [aq,a*,] — Sq,q>. The 
Hilbert space of the total system consisting of small system and field, is ,^^ = ,J>^ (8) J^, and we simply write Hg 
and H^ for the operators Hs ® 1 and 1 ® H^ acting on J^^. The coupling between field and the small system is 
assumed to be linear in the creation and annihilation operators and it can hence be written in the form 



Hi, :=. ^ A ® Hcj^f), ^4) := J2 {'^HqK + '/'f (9)S) (1-2) 

iei qeA' 

where Di = D* are self-adjoint operators on Jfs and (fyf- are functions (form factors) to be specified. X is a finite 
index set. 

The total Hamiltonian of the system is hence, with a coupling strength A e R, 

H^ := Hs + H^ + XH,i,, on JT^. (1.3) 

A standard application of the Kato-Rellich theorem states that, if (cj'^)^^/^(/)f^ e P{A*), then H^ is self-adjoint on 
the domain of H^. 

Initially, the field is in a Gaussian state that will be called the 'reference state': The density matrix pp ' of this 
reference state is specified by the correlation functions of smeared field operators, i.e. by 



Tr, 



ef.A, 



p^"^'^$(^i)...$(^„) , (1.4) 



where ipi are functions on A* and Tr^ is the trace on J^^. 

Our main assumptions on this reference state p^ ' are that it is 



1) stationary w.r.t. the decoupled dynamics, i.e. 

2) gauge-invariant, i.e. all correlation functions l|1.4l l that involve an odd number of field operators, are zero. 

3) Gaussian (also called "quasifree"), i.e. the higher correlation functions are related to the two-particle correla- 
tion function via the Gaussian relation 



Tr,, 



pairings IT (i,j)e7r 



where a 'pairing' tt is a set of pairs (i, j) with the convention i < j. 

By standard theory of Gaussian (or quasi-free) states, the above properties imply that the state p^° ' is completely 
determined by a positive density function < ri^{q) < cx) via the relations 



Trp 



pf'^<P{^Mv') = (^', 77^^) + (^, (1 + 77^)^') (1.7) 



cf,A 



The invariance of pp ' under the free field dynamics is ensured by the commutation relation [uj , rj ] ~ 0. Since 
the field is a finite collection of harmonic oscillators, the reference state p^^/ is a well-defined traceclass density 
matrix. The material in this section is completely standard and we refer the reader to e.g. |7, 4] for details that 
were omitted here (note however that these texts deal with infinite volume A from the start and hence they are 
necessarily more involved technically). 

1.3 Thermodynamic limit 

As long as A is finite, we cannot expect the system to have good ergodic properties and hence we will perform 
the thermodynamic limit as a first step. By the thermod3n-iamic limit A /* Z"^, we mean that L /• oo, hence the 
volume A tends to Z"* and the set A* tends to T"^. As we will see, the influence of the reservoir on the dynamics of 
the subsystem can be expressed entirely in terms of the correlation functions 

/^(t) := TrRi/,^°f''^$(e'*-"0f)$(0A)] (l.g) 



yr^ ' •! ^ •Til' \-ri IV-" ' '/ )^ vj 



^ r;^e'*- 0f ) + {<l>t (1 + ^^)e-"" <) (1.9) 



Note that j^M) = fji{—t) by stationarity of the state pp ' . To discuss the thermodjrnamic limit of the small 
system behaviour, it suffices to ask that the /,A (i) converge 

A,(t)= lim /5(i), (1.10) 

uniformly in t on any compact set, and that sup^ \fi.j{t)\ < oo. We recall that a density matrix on J^s is a positive 
traceclass operator, i.e., belonging to ^i(J^), whose trace is 1 (of course, since J^s is finite-dimensional, any 
operator is traceclass). In what follows, we let Tir stand for the partial trace over the field degrees of freedom, 
mapping density matrices on M'^ into density matrices on Ms ■ 

Lemma 1.1. Assume that the fi^j{t) converge to hounded functions fi.j{t), uniformly on compacts (i.e. \\.1Q) ). Then, the 
thermodynamic limit 



ps,t ■= lim TrR 

A/'Z'i 



e-««. 



{ps.,®pt^)e^'"^] (1.11) 



exists for any initial density matrix psfi on ,^. 

The proof of Lemma [TT] is given in Section l3.1.2[ 



1.4 Markov approximation 

The Markov approximation to the model introduced above amounts to replacing the correlation function fij {t) by 
a multiple of S{t) (no memory). It can be justified in the weak coupling scaling limit A — >^ 0, if one rescales time as 
t — > X^^t. We state this important result precisely in Section l331 For now, we just introduce the precise form of the 
Markov approximation since one of our assumptions refers to it. First, we introduce the left and right multiplication 
operators MLe(A), M^i^A); 

Mi^e{A)S -.^ AS, Mjii{A)S -.^SA*, A,Se^{^s) (1-12) 

Then we set 

,-^ /"^ , s , ■. { fii{t) if ki = Le 
L:= ^ / dte''^'^^"^^Mk,iiDj)e''''"^^"^^MkAiD^)l (1.13) 

where ad(iJs) = [Hs, ■] and the integral over t is well-defined by integrability of fi,j{), which will be assumed 
below. Finally 

1 pt 
L-~ lim - / dse''"''^(^^)Le-'"^^(-f^^) (1.14) 

t^OO t Jq 

where the limit t -^ oo exists since Hs has discrete spectrum. Note also that L commutes with ad(iJs) as follows 
from the spectral averaging in l|1.14l l. As is discussed in many places, the Lindblad operator L generates a contract- 
ive semigroup e*^^, t > on 3ii{Jfs) that is trace-preserving and positivity preserving. In other words, e*^ maps 
the set of density matrices on ,J^ into itself. In the above formulas, we denote time by the gothic symbol t to em- 
phasize that it corresponds physically to a rescaled time. Indeed, the Lindblad operator L describes the dynamics 
on long time scales, see Section |33l Lindblad operators were first introduced in Il22l , an excellent exposition on the 
properties of L and its derivation from microscopic models can be found in [20l . 

1.5 Result 

We need an assumption on the decay of temporal correlations of the 'free reservoir correlation functions'. 
Assumption 1.2 (Decay of correlations). Recall the correlation functions fi,j introduced in (|1.10l l. We assume that 

/ dth{t)<oo, where /i(i) := ^ ||A|||P,|||Aj(t)| (1.15) 

The second assumption concerns the Lindblad generator L, defined in Section [l~4l 

Assumption 1.3 (Fermi Golden Rule). The operator L has a simple eigenvalue at 0. All other eigenvalues lie in the 
region {z e C | Re z < — gap^} for some gap^ > 0. 

Obviously, Assumption ll .3l and the fact that e*^ preserves density matrices, imply that there is a unique density 
matrix, pg , such that Lp^ — and 

||e*^ - \p^){l\ II < CiC-s^P^*, for all t > and some Cl < oo (1.16) 

where || • || is the operator norm of operators acting on ^i(J^) and we use the notation \A) {A'\ to denote the rank-1 
operator that acts as 5 -> {Tt[{A')*S])A with S,A e ,^i{,ytfs) and A' e ^(J^). For us it is more convenient to 
define a characteristic time ii > 1/gap^ such that 

||c*^-|pi')(l|||<c-'/*-, fort>ii (1.17) 



Conditions that imply Assumption !! .3l have been discussed extensively, see e.g. [;T2ll28|. Here, we prefer to give 
a (rather generic) example where the Assumption ll.3l can be checked very explicitly: Assume that the Hamiltonian 
Hs is non-degenerate, hence its spectral projections, -P(e), e G spiJg, are one-dimensional. Then Assumption [L3] is 
satisfied if and only if the continuous-time Markov procesqj with (finite) state space spiJg and jump rates 

rate{e ^ e') = V /,;,,(e - e') Tr[P(e')i?*P(e)A], where f,,,{u) ^^ f dte^'''^ f,,,{t) (1.18) 

is ergodic. This in turn can be checked by the Perron-Frobenius theorem: a sufficient condition for ergodicity is that 
for any two eigenvalues e, e', there is a path eg, ei, ..., e„ with Go = e, e„ = e' such that, for alii, rate (e^ -^ Si+i) 7^ 0. 
We are now ready to state our main result 

Theorem 1.4. Assume that Assumption \1.2\ and Assumption \1.3\ hold and let ps,t bs defined as in Lemma \T?I\ Then, there is 
a Xq > such that for all A satisfying < |A| < Aq, H'e have 

lim ps,t = pT (1-19) 

where the "invariant density matrix" pg"" = Pg"^(A) does not depend on the initial state ps,o- Moreover, pg"" is a small 
perturbation of pg, the invariant density matrix predicted by the Markov approximation; 

WPs"" - PsW -^ ^ fls A^O (1.20) 

In (|1.20l l, II • II is the operator norm of operators acting on ^i(J^) (although it does not matter since J^s is 
finite-dimensional) . 

To quantify the speed of convergence towards the steady state Pg"", we need to know the decay properties of 
the function h{-) that was introduced in Assumption ll.2l Let ({■) be a nondecreasing function on M+ satisfying the 
conditions 

l<C{t + t') < C(i)C(i')> for any i, t' e M+ (1.21) 

We assume that this function governs the decay of the bath correlation function h, in the sense that 



dth{t)Cit) < 00 (1.22) 

/o 

The case where (^{t) can be chosen to be exponentially increasing, is particularly simple but introduces a complica- 
tion to the statement of the following result. Therefore we exclude this case explicitly by demanding 

dte''^X{t) < 00, for any k > (1.23) 

Proposition 1.5. Assume the conditions of Theorem \1.4\ Let Q be a non-decreasing function as above, satisfying (12. 221 
\1.22Vil.23i , and let i^ be chosen such that HI. 17b holds. Then, for \X\ small enough, 

\\ps,t - PrW < exp (-^7^^^) + «(l^r) (C (1^)) > for any t > i^X-' (1.24) 

Note that Proposition 11.51 makes no claim about the reduced dynamics ps.t for short times t < A^^t^. The 
restriction to long times is natural since, for times shorter than A^'^tL, the exponential decay of the semigroup is 
not yet visible. For those times, ps.t is however well-described by the Markov approximation, see Theorem l3.5l On 
the RHS of l|1.24|l . the time t appears essentially in the combination X^t/iL. As far as the first term is concerned, 
this is natural since that term originates from the Markov approximation, i.e. the temporal decay embodied in that 
term takes place on the macroscopic time scale ^ A^^l^. The second term, however, comes from the slow decay of 
the reservoir correlation function h{t) on the microscopic time scale, and as such it is not clear why that decay gets 
prolonged to the macroscopic scale in l|1.24|l . The estimate in that second term is indeed far from optimal (note also 
the weird factor '2' multiplying t^) and this is due to the generality of our result. If, for example, one assumes that 
C{t) ^ |t|", a > 0, then one can state a sharper and more explicit bound. 



^Since L commutes with ad(iJs) and preserves positive density matrices, it sends the set of density matrices 
diagonal in iJs -basis into itself. Since these diagonal density matrices can be identified with probability measures 
on spHs, e'^ determines a Markov process on spHs, namely the one defined by the rates l|1.18t . 



1.6 Discussion and comparison with earlier results 

1.6.1 Restriction to confined systems 

Our result is suited for confined small systems. We explain this in more detail and we distinguish essential as- 
sumptions from those made for convenience. 

A: The assumption that the 'atom' Hilbert space J^s is finite-dimensional, seems not crucial to us. Atoms 
with an infinite number of energy levels (like the harmonic oscillator) should be treatable with the same 
technique. A complication that does arise in such infinite-dimensional atoms is that the relaxation of the 
Markov semigroup is in general not exponential since, in the absence of very energetic field quanta, the 
atom needs a large time to cascade from a very energetic level to the low-lying levels. We believe however 
that this can be remedied by a change of norm on (a subspace of) J^ that renders the relaxation exponential, 
at least for a certain class of interaction Hamiltonians. 

B: The restriction to atom-bath couplings that are linear in the field operators is for notational simplicity only. 
One can study quadratic coupling in the same way. Coupling terms of higher order do not yield a well- 
defined Hamiltonian for bosonic baths, although they are well-defined for fermionic systems. In that case 
(fermionic baths with coupling of order at least 3) one has to use sign cancellations to control the Dyson 
expansion (this is done e.g. in 1161 ) and in such cases an operator-theoretic treatment might be favorable. 

C: The real assumption that excludes application of our result to extended systems is Assumption ll.2l and more 
concretely, the sum over i^j el. For an extended system, the simplest translation invariant coupling would 
be of the form 

i/int = / dg¥>(g)e'«^a,+^((7)e-''^a* (1.25) 

= ^ |x)(a;|®$((p,), with^,(g)=e'«Xq) (1.26) 

where we have taken k = U^ . The expression on the second line is of the form (|1.2l l with the index set 
1 = 1/^. Even though one could demand that the correlation functions are integrable in time in the sense 
that 

r+oo 

sup/ dt/la;,x'(i) <00, (1.27) 

x^x' JO 

with 

icf,A I / ito;. 



h^^^,{t) := lim TrR p™''^$(e'*"^,)$(^,0 



(1.28) 



then still Assumption 11.21 cannot hold because of the sum over x,x' E JJ^. In fact, the appearance of the 
double sum is artificial and one can arrange to have a single sum, and moreover, h depends on the difference 
x' — X only. Hence, Assumption ll.2l would boil down to 

/>CXD 

/ di^VxW<«3 (1.29) 

Jo 

and this assumption cannot be satisfied for any interaction Hamiltonian of the type (|1.251 . 

1.6.2 Interacting reservoirs 

Models where the heat bath is not free, i.e. it is made up of a genuinely interacting system, are a far dream at this 
moment. However, we would like to draw attention to the fact that, in contrast to earlier results, our method does 
not exclude such reservoirs per se. Indeed, the important ingredient of our analysis is a temporal decay condition 
on the reservoir correlation function. This condition is stated in Assumption l|2.2l l, and, for free reservoirs, it is 
satisfied provided that Assumption [L2] holds. The huge challenge is of course to prove such a condition for an 
interacting system. First steps in this direction have recently been taken in L24J. 



1.6.3 Algebraic quantum dynamical systems 

In the literature on the subject, mixing properties are mostly investigated in a more general framework, allowing 
for initial states that are not factorized (but still local perturbations of pj^^) and treating observables that depend 
on the field as well (since Theorem 1 1 .41 deals with the reduced dynamics, we get information on observables of the 
small system only). In particular, one usually studies the system in the framework of C* or M^*-algebras, in which 
the concepts "ergodicity" and "mixing" have a natural meaning, inherited from the theory of d5mamical systems. 
For an introduction to these matters, we refer to [T'4\. It is straightforward to extend our approach such as to prove 
mixing in the above sense, but since this asks for more notation in Section|3l we have opted not to do so. The same 
remark applies to the study of multitime-correlation functions of small system observables. Our technique shows 
that these correlation functions are perturbatively close to correlation functions calculated within the Markovian 
approximatiorQ see also [llj. A draw^back of our technique with respect to the algebraic approach is that, in the 
case where p^^ is a Gibbs state, it is not immediately clear that the invariant state p™" is the restriction of the 
coupled Gibbs state to the small system. However, if one extends the class of initial states as suggested above, this 
does immediately follow. 

1.6.4 Comparison with earlier results 

One should distinguish between the case where the Gaussian reference state of the field has a non-zero density 
(temperature) in the thermodynamic limit, i.e. liui/^y^^d |A|^^ X^qga* V^il) > 0' ^^ ^'^^- ^^ ^^^ latter case, the field 
is essentially in the vacuum state and the approach to a steady state is related to the question whether the ground 
state of the coupled system (assuming that it exists) is the only bound state and whether the rest of the spectrum is 
absolutely continuous. These questions have been extensively studied in IIT5l [Tll3l[T4l. In one sense, our results are 
sharper than those quoted: they cover cases where the coupled system has no ground state, yet there is approach 
to a steady state for the small system. We do not explain nor develop these issues further here, but rather postpone 
them to a subsequent paper. However, the quoted results are stronger in the sense that they allow for the confined 
system to have continuous spectrum above a ionization threshold. 

If the field has a positive density, the prime example is of course the case where the field is in a thermal state 
at non-zero temperature, then the only results that we are aware of, rely on complex deformations of Liouvillians. 
One either uses complex translations or dilations. To streamline the discussion, we note that one can rewrite the 
correlation functions fi.j{t) as 



Ujit)^ / dwe'*"/,,,(w), /,,jH:=(v?,,e//H,Vj,e//(w)),i^ (1.30) 

such that fi.eff, the effective form factors (in the thermodynamic limit), are functions from M to some Hilbert 
space S^ that emerge naturally if one follows the operator-theoretic approach to the problem. They are often called 
"effective form factors"(effective because they incorporate the density function ( of the reservoir). In the physical 
literature on the subject, the function fij (w) is often called the 'spectral function'. 
The first result on RTE, due to flTllTHl , proceeds by assuming that 

• The function uj n- fi,j{oj) is analytic in a strip of width 70 such that uj h^ fi.j{^ + 17) is in L-'^(R, doj) for 
< 7 < 70. 

This of course corresponds to exponential decay of ,fij{t). This result has been improved in HISl where 
analyticity is replaced by demanding that /^ j{-) is in C^, implying fi j{t) ^ |t|^^. A related approach is found in 

The approach via dilation analyticity has been pioneered by [2j. There one assumes that 

• the function uj 1-^ /^ ., (e'('''S""''''w) is in L^{R, duj) for < 7 < 70 (this is dilation analyticity) 

• /»j(e'("'en<^)''w) < |a;|i+" for some a > 0. 



^Yet, they are qualitatively different, since in the Markovian model, correlations decay exponentially, whereas 
at finite A, the speed of decay is in general not faster than the decay of the correlation functions fi,j (t). 



By deforming the integration contour M in l ll.30|l into c '''M- U c'''^M+ , one realizes that this implies that 

\Aj{t)\ < const t-(2+") (log f)™-*' (1.31) 

and hence this case is covered by our result. 

1.7 Strategy of the proof 

Our proof is based on a polymer expansion in real time. In the context of classical stochastic dynamics, such 
expansions were successfully applied in e.g. Il5l l25ll , and in the case of classical deterministic dynamics in L19J . For 
the case at hand, a similar strategy was pursued in l27l . In the following Sections [l.7.1[|1.7.2l and [m3l we introduce 
the rough ideas. 

1.7.1 Markovian approximation and leading dynamics 

We discretize time t — Nv where i^ is a macroscopic time unit v = X^^l, with I a A-rndependent number that 
could actually be chosen £ ^ 1. Then, we write ps.t = Z^ps.o where ps,t is the reduced time-evolved density 
matrix and Z^ could be called the 'reduced evolution operator'. The idea is that T = Zm=i can be analyzed quite 
well, at least for sufficiently small coupling A, because in that regime the Markovian approximation (Section II. 4t 

can be justified. Indeed, we will state in Section l33l that T is well-approximated by e^ '^^ , with L the Lindblad 
generator (also mentioned in Section [L4l l. This is not proven in the present paper since the proof is well-known 
in the literature. For now, we view T as the leading dynamics. An important consequence of the fact that T is 
close to 6^ ^^ and of Assumption 11.31 is that we can establish that the operator T has a simple eigenvalue 1 (this 
eigenvector is the 'steady state' pa^) and the rest of the spectrum lies in a circle with radius 1 — <? < 1. Since T is trace 
conserving, Tr Tpgo — Tr ps.O/ the 'right' eigenvector corresponding to the eigenvector 1 is the identity 1 e ^(^ ), 
hence we have the spectral decomposition 

T^R+{l-R)T, R=\pI){1\, ||(l-i?)T"||<C(l-g)" (1.32) 



This property of T is proven in Section 13.31 by simple perturbation theory (with c'^ ''^ being the 'unperturbed 
object') , but it is introduced already in Section |2] as an assumption. 

1.7.2 Polymer representation 

If the reduced dynamics Zm were exactly Markovian, we would have Zjv = T^ , i.e. Zjv could be called a 'quantum 
Markov chain'. However, this is of course not the case and as N grows the difference between Zm and T^ becomes 
important. We represent the corrections to T^ by 'nonmarkovian excitations' that are localized in time. For ex- 
ample, 

Z2=r2 + r[E^(B(l,2))] (1.33) 

where W{B{1, 2)) is an operator on ^i(^) ® ^i(^) that should be thought of as localized in the macroscopic 
times 1 and 2 (actually, in macroscopic time intervals [0, i^] and [v,2v]). The operation T[-] is a time-ordering; it 
converts W{B{1, 2)) to an operator on ,^i(J^), such that it is on the same footing as T (see the full definition in 
Section IZ21 We are actually abusing the correct definition slightly in the present section). For Z-^, we get 

Z3 = T^ + Tr[E=(B(i, 2))] + rwiB{2, 3))]r + r[E=(s({i, 2, 3})] + r [ie"(b(i, 3))r(2)] (1.34) 

where, for a general set of macroscopic times A, W{B{A)) denotes the (irreducible) excitation that is localized in 
the elements of A (it acts on the |A|-fold tensor power of 3§i{Ms))- Since, in the rightmost term, the excitation 
W{B{1, 3)) is localized in times 1 and 3, and the T-operator represents the leading dynamics in the second time 
interval (hence the '2' in T(2)), we need to squeeze T{2) in between the excitations at times 1 and 3. 



For general N, the resulting expression for Z^ is 



AeVo\(N) 



V rG/jv\Supp^ 




W{B{A)) 



AeA 



(1.35) 



where the polymer set Pol(A^) is the set of nonempty collections A of disjoint subsets Aoi In = {1,2,..., N}. To 
analyze this polymer expression, we use two tools: bounds on the excitation operators W{B{A)) and a Feynman 
rule. 

Bounds We will bound each term in the sum lll.35l l in operator norm by n^eyt ll^'^(^(^))ll# ^vhere the norm || • ||# 
is defined in Section r2.2.2l The r(r)-operators do not show up in these bounds since they have norm 1. We will 
require that ||E'^(i3(A))||^ ^ el'^l for some small parameter e and, moreover, that ||E'^(_B(yl))||^ decreases as the 
macroscopic times, i.e. the element of A; are further apart. This decrease as a function of temporal distance is a 
consequence of Assumption ll.2l but in Section|2]it is introduced as Assumption l2.2l In Section l3.1.1i we prove how 
Assumption ll.2l implies Assumption l2.2l 

Feynman rule It is not hard to see that the bounds given above, when summed over the different terms in (|1.35ll 
lead to a too pessimistic bound on Zm . Even if we restrict to sets A whose elements are consecutive integers (which 
is essentially justified because of the temporal decay), then we still get an exponentially diverging bound, of order 
e'^'^^, for some constant C > 0. To improve our bounds, we use a Feynman rule (one could also call it a Ward 
identity) that is a consequence of conservation of probability of the dynamics Z^, to be explained in Section l23l 
In our general polymer expansion, this Feynman rule implies that, for every uninterrupted string of T() operators 
that follows a set A, we can insert the spectral projection (1 — i?) in front of the string of T's. By l ll.32l l, this yields 
exponential decay in the length of the string. This is illustrated in Figure[T](the sets Hook(y4) will be defined later). 
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Figure 1: An example oi a A with A = {Ai,A2,Aj,}. In the picture, A^ = 20, and Ai == {3,4, 6}, A2 = 
{5, 10, 11, 13}, A3 = {16, 17}. The exponential decay is on the string of times that are covered by the dotted lines. 
These strings are Hook(yli) = {7, 8, 9}, Hook(A2) = {14, 15}, Hook(yl3) = {18, 19, 20}. These are exactly the times 
between max A for some A and the next-in-time element of some other set A' . 



Armed with the Feynman rule and the bounds on W^, we can now perform the sum over all terms on the RITS 
of l ll.35|l , resulting in 

||Zw-r^||=0(e) ^ \\ZN-R\\^C{l-g)'' + 0{e) (1.36) 



By very similar arguments, one can see that Z^, for large N is e-close to a one-dimensional projector. Indeed, by 
the exponential decay following any excitation, all contributions to Z^r — T^ tend to be localized at times close to 
the final time N. This means that they start with a string of T's of length 0{N); such a string is equal to R, up to an 
error of order (1 — g)'^'^-'. Hence, up to a vanishing error, all contributions get multiplied by R on the right, and 
consequently they are of the form \S) (1| for some S € ^{J^s)- This means that also the Hmit limAr/<oo Zn is of this 
form. By conservation of trace and positivity it then follows that 



lira Zn = 

N/'oo 



pDiM 



for some density matrix pg"" : ||pg" 



pI\ 



= 0{e), 



(1.37) 



These arguments are presented in Section l2!6l 



1.7.3 Dyson expansion 

The Dyson expansion is introduced to prove the bounds on j|E'^(i3(A))||^ discussed above. This is done in Section 
13.21 It is also the standard tool to prove the weak coupling limit, Theorem l3.5l 

As we will do later on in the proof, we assume for simplicity that the set X has just one element, such that 
we can drop the index z S I and simply write H^^^ = D ® $((/3*). For any operator O on J^f^, let 0{t) = 
Qit(Hs+Hj^)Q^-it{Hs+Hj^) ^ and we abbreviate $(f) = (<l>((/)^))(t), then we can write the Duhamel expansion (the 
convergence of the series is easily established) 



gitad(ffs) 



ps,t ■■= lim e'*^^TrR[e-*^^(ps.o®PR''^)e"«^1e-"^^ (1.38) 

= 1%, E (-iA)"^"(iA)"- j dsi...d5„,, j ds;...d<^^ 

= J2 (-iA)"-(iA)"-^lim^TrR [^s,,,J . . .^s,Ms,)p'^''^^s[M4) . . .^s'^^ 



ds J ds' D{s,,,J...D{s2)D{si)ps,oD{s[)D{s'2)...D{s'^J 

0<si<...<5„, <t 0<s'<...<s' <t 

Note that the operators on the last line act trivially on J^ji, and hence they are independent of the volume 
A. The expression on the one-but-last line is a riLe + ?T^Ri-point correlation function corresponding to the Gaussian 
(quasi-free) state p'^ ' . Hence, by the Wick theorem, we can expand this correlation function into sums of pairings 
of the nLc + njii of products of two-point correlation functions /^(•), as in formula jl.bi . Each term in the sum is 
determined by a pairing of the ulo + nRj-times, it is called a "diagram". 

Starting from the Dyson expansion, we first identify which terms in that expansion make up the "excitation 
operators" K'^{B{A)). This is particularly intuitive. For example, the operator E''{B{{t, t'})) is built by all terms 
(diagrams) in the Dyson expansion whose times s, s' fall into the domain (/[t — 1, r] U i/[t' — 1, r'] and such that at 
least one pair in the pairing connects the two intervals, i.e. it has one of its time-coordinates in each interval. For a 
general set of times A, the operator W{B{A)) is made up by diagrams such that the set A is connected by the pairs 
in that diagram Two examples will be given in Figured 

2 Polymer model 

In this section, we start from a discrete-time dynamical system and we derive the approach to a steady state, given 
some assumptions that will be justified in Section|3l Apart from the first paragraphs, the discussion in this Section 
is independent of the setup given in Section [l] In particular, it could be applied without any change to other 
models, hinted at in points A and B of Section [l.6.1l 

2.1 Reduced dynamics and excitations 

We define the propagator U{t), implementing the dynamics between macroscopic times r — 1 and r, with r G N, 
and acting on joint density matrices psR/ 

^(•.^) .^ gi,.rad(ff^)g-i..ad(H^)g-i,.(r-l)ad(i/^) (2.1) 

where we have chosen the macroscopic times to be related to the microscopic times by a scaling factor v that will 
be fixed in Section |3l depending on details of the model. The total dynamics (in the interaction picture) up to 
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microscopic time N is then given by 

UiN) . . . U{2)U{1) = ei''^ad(ffi^)g-iz.Arad(//i^) (2.2) 

Furthermore, we define a projection operator P; 

PpsR-=iTrRPsii)(^Pn''' (2.3) 

A distinguished role is played by the reduced dynamics T{t), acting on density matrices ps on J^. It is defined 

by 

iT{r)ps) ® Pt^ = PU{r){ps ® p'^''^) (2.4) 

Note that T{t) is independent of r by the stationarity property lll.SI l, but we still write the dependence on r for 
bookkeeping reasons that will become clear in Section lZ2l It is convenient to abuse the notation and let T{t) stand 
for T(r) (g) 1 as well, such that it acts on joint density matrices. (This abuse of notation appears only in the present 
section.) Note that, as operators on joint density matrices, 

T(r)(l-P) = (l-P)T(r)^0 (2.5) 

Define the 'excitation operator' 

P(r) := U{t) ~ T{t) (2.6) 

and note that 

PB{t)P = (2.7) 

Ultimately, we are interested in the reduced d3n-iamics: 



^N 



PU{N) . . . U{2)U{l)P = Pe~i'^N^d{Hi)p (2.8) 



We will next insert the decomposition U{t) — T{t) + B{t). We note that the definitions of U{t),B{t) are given 
only in finite volume (|A| < oo). However, in Section |Z41 we will obtain some concepts (including Zn) that do 
admit a thermodynamical limit. 

2.2 Correlation functions of excitations 

For notational purposes, it makes sense to define the following: Let V be an operator in (^reiN^r where In = 

{1, . . . , N} and each ^T- isacopy of ^ = ^i(J^). Define the 'time-ordering" T as a linear operator (g)T-g/„ ^(.^7.) — >■ 
^{^) : V I-)- T[V] as follows. For elementary tensors V ~ (SDre/jv^r where Vt G ."^t, we simply put 

T[V]:^Vn...V2Vi (2.9) 

and then we extend T by linearity to the whole of ®i-e/jv^(-^r)- Now, take a subset A c In and define the 
operator 

E(B(A)) : ^A -^ ^A, where Ma := «)reA^r (2.10) 

as follows. Let A ~ {ti, . . . , r™}, with the convention that ti < T2 < . . . < t„j, choose operators Sr, S'^ in Mt.S^J^, 
respectively (here, S^.'^ is the dual space to ^r) and \e-iSA — ®T^AST and 5^ = ^tgaS'^ be elements of Si. a and ^^ 
, respectively . We define the operator E(P(A)) by giving its 'matrix elements', namely 

{S'A.nB{A))SA) :- Tr[5;^^P(T„0(|5,,J(5;_J®lR)... (2.11) 

. . . P(r3) (|5.3)(5;j ® Ir) P(r2) (|5.,)(^;| ® 1r) P(ri)(5., ® p^^^'^)] (2.12) 

where (•, •) on the LHS is the pairing between S/?\ and ,^^ and, on the RHS, the notation |A)(A'| for operators 
on i^ ^ ^i(J^) was introduced following Assumption 11.31 Note that E(P(y4)) = whenever A has only one 
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of, A 



element, since the operator PB{t)P vanishes on tensor products of the form S^p'^ ' . Furthermore, the correlation 
functions corresponding to sets A and A + t are copies of each other, but acting on different spaces {^a vs. Ma+t)- 
This follows from the stationarity of the reference states pp ' under the free reservoir dynamics. By expanding 
U{t) — T{t) + B{t) for every r in the expression for the reduced dynamics (|2.8b , we arrive at 



■^N 



-Er 



Ac/a 



T{r) 



miB{A)) 



,re/jv\A 



(2.13) 



On the RHS, the first tensor acts on Mjj^\a and the second on ^a- The time-ordering T makes sure the operators 
are 'contracted' in the right way. Note that the order in which we write the tensors in expressions like ||2.13|| does 
not have any significance. Instead, the 'legs' of the tensor product on which the operators act are indicated by the 
arguments t and A. 

2.2.1 Connected correlation functions 

The "connected correlation functions", denoted by W{B{A)), are defined to be operators on ^a satisfying 



i{BiA')) = Yl ^'iB{A)) 

A partitions of A' \AeA j 



(2.14) 



The tensor product in this formula makes sense since S^-a' = ®Ai^a^^a whenever ^ is a partition of A'. Note that 
this definition of connected correlation functions reduces to the usual probabilistic definition when all operators 
that appear are numbers and the tensor product can be replaced by multiplication. Just as in the probabilistic case, 
the relations (|2.14) l for all sets A! fix the operators W{B{A)) uniquely since the formula (|2.14) l can be inverted. 
With this machinery in place, we can write a neat expression for the reduced dynamics Zj^; 



■^N 



T 



N 



A&Vo\(N) 



with Zn{A):^T 



T{t) 



re/jv\SuppA 




W{B{A)) 



AeA 



(2.15) 



where the polymer set Pol( A^) is the set of non-empty collections A of disjoint subsets A of In, and Supp.A — yjAeAA. 
Formula | |2.15l l follows from | |2.13| | by substituting l|2.14b since, obviously, any polymer A\s a partition of Supp^. 
The term T^ in l|2.15l l originates from A = in l l2.13l l. 



2.2.2 Norms on ^{Ma) 

We introduce a norm on the spaces S§{i^a) that is appropriate for our bounds. Any operator E on ^^ can be 
written (usually in a non-unique way) as a finite sum of elementary tensors 



E = Y.E,, 



(2.16) 



We define 



i^ii#^=/^^5jE 11^-^1 



(2.17) 



where the infimum ranges over all such elementary tensor-representations of E. The sole purpose of the norm 
II • ||# resides in the following estimate: for any A E Pol(iV) and any operators Cr on ^r, ^ve have 



r 



E^B{A)) 



Cr 



\AeA / \Te/iv\SuppA 

This follows immediately from the definition (|2.17b . 



AeA re/iv\SuppA 



(2.18) 
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2.3 Unitarity 

We now examine the consequences of the conservation of probability, i.e. unitarity of the propagators U{t). 

TrB(r)psR = 0, for any psR e ^i(^s ® ^^r ) (2.19) 

This follows from the fact that 

Tr B{t)psr = Tr U{t)psr - Tr(r(T) l)psR = Tr U{t)psr - Trg TrR U{t)psr = (2.20) 

Let R be a one-dimensional projector on ^ of the form R = |ps)(l|/ for some density matrix ps e .^i(J^) (i.e. 
satisfying Tr ps = !)• It follows that 

i^^{^A^^..A} ® i?(max A)) E{B{A)) = 0, (1,^(5?^,_^) ® i?(max A)) E=(B(A)) = (2.21) 

where R{t) is a copy of R acting on ^t (here used for r = max A) 

2.4 Infinite-volume setup 

As remarked already, the above setup makes sense in finite volume A only. However, the operators Z^ , T and 
the "correlation functions" M{B{A)), E''{B{A)) have well-defined thermodynamic limits, since these operators act 
on (tensor products of) the system space only. More importantly, the relations l|2.13t , | |2.15|I and l|2.21|l will remain 
valid in the thermodynamic limit. For the sake of explicitness, we put this in a lemma 



Lemma 2.1. The correlation functions E{B{A)),W{B{A)), as introduced in Section |Z2] depend on the volume A via the 
reference state p^ ' , and hence we should denote them by E^{B{A)),E,^''^{B{A)). However, the infinite-volume correlation- 
functions 

E{B{A)) := lim E^{B{A)), EHBlA)) := lim E"'^(BM)) (2.22) 

exist and satisfy the equality l|2.21|l . 

In the remainder of Section|2l when writing 'E{B{A) ) , E'^ {B{A) ) we always mean the infinite-volume quantities. 
The proof of Lemma IZT] is analogous to the proof of Lemma Fl. II which is contained in Section l3.1.2l 

2.5 Abstract result 

We now state two assumptions that allow us to prove the convergence to a stationary state. In Section 3 we verify 
these assumptions for the model considered in Section 1. 

Our first assumption concerns the decay of truncated correlation functions. We demand that A: -point trun- 
cated correlation functions have a sort of tree-graph decay reminiscent of high temperature Gibbs states. Let 

A = (ti , . . . , Tfe) with n < Ti+i and put 

fe-i 
distc(A) = distc(T) := H C(|r,;+i - n\) (2.23) 

with ( a nondecreasing function satisfying i ll. 2111 and l| 1.231 

Assumption 2.2 (Summable correlations). There is a constant e < cjo such that 

sup Y. e-\^\distc{A)\\WiB{A))y<l (2.24) 

max A = TO 

uniformly in N, 
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The next assumption expresses that the operator T has a well-defined leading order contribution. 

Assumption 2.3 (Dissipativity). The operator T has a simple eigenvalue 1 corresponding to the one-dimensional 
spectral projector R — jpg ) (1| with pE^ a density matrix. Moreover, the rest of the spectrum lies inside a disk with 
radius 1 — g with .g > 0. 

This assumption obviously implies that, for some constant Cg < oo, 

||T"(1 - i?)|| < Cg(l - .9)", for anyn e N (2.25) 

Our main theorem states that under these assumptions, the system state approaches a unique limit as A^ — >^ oo 
Theorem 2.4. If Assumptions\2.2\and \2.3\ hold with e sufficiently small, then there exists a density matrix p™" such that 

lim ZA, = |pr>(l| (2.26) 

where 

\\pl-pr\\^0{e), e^O (2.27) 

The speed of convergence is estimated as 

\\ZM~\pr){M\\<Cg{l-9r + 0{e)^y e^O (2.28) 

2.6 Proof of Theorem D 

2.6.1 Summation of the polymer series 

We start from the representation 

Zjv=r^+ Y. ^^(-^) (2.29) 

AeFol{N) 

and we use the conservation of probability (Section |Z5t with 

R=\pI){1\, (2.30) 

where pg is the unique invariant state of the map T. It follows that 



Zn{A)=T 



(g)E^(i3(A)) (g) (g) T(r)(l-i?) (g) (g) T{r) 

\AeA I \TeHook(^) / \/jv\(Supp^UHook(.A)) 



(2.31) 



where the set of times Hook(yl) is determined as follows. For each set A E A, pick the smallest time in Supp^l \ A 
that is larger than max^ and call this time THook(^)- If there is no such time (which happens for exactly one A, 
namely the one for which max A = max Suppyl), then set rHook(^) = A^ + 1. We also define the sets 

Hook(A) := {maxA+l,maxA + 2,...,THook(A) - 1}, Hook(yt) := |J Hook(yl) (2.32) 

AeA 

(in case max A + 1 = T^ookiA), we set Hook(A) :— 0.) In words, the set Hook(74) is simply the set of times between 
the latest time of A and the next element that belongs to some A' e A. See Figure [l] for a graphical representation 
of a term Zn{A). 

At this point, there is no evident relation between operators Zn{A) for different N. However, let m{A) = 
min Supp(^) denote the earliest time included in the polymer A and take a A G Pol(A^') with mlA) > N' — N, for 
N' > TV. Then, the polymer A = {Ai ,... , A^} can be converted into a polymer in Pol(7V), namely A — {N' — N) := 
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{Ai-{N' -N),...,Ara- [N' ~ N)} where A, - [N' ~ N) ■- {r - (TV' - iV), r e A,}. By putting the projector R 
on the right of Zn{A), we obtain then that 

Zn'{A)R = Zn{A-{N' -N)})R, iim{A)>N'-N (2.33) 

This is a straightforward consequence of the translation invariance of the correlation functions. Hence the differ- 
ence between Z^R and Z^'R is made by polymers that are "too extended" to be fitted into In {m{A) is too small); 
we have 

{Zn' - Zn) - (T^' - T^) (2.34) 

J2 Zn'{A)R+ Y. ZN'{A)il-R)- Yl ZNiA){l-R) 

; Pol(Ar') 

<N' - N 

Lemma 2.5. If e is small enough, then 



A G Pol(Ar') A G Pol(Ar') A G Pol(Af) 

m{A) <N' - N 



Y \\ZNiA)\\=0{e), 6^0 (2.35) 

AePoI{N) 

and 

lim y ||Zjv(-4)(l-i?)|| -0 (2.36) 

.AePoi(JV) 

Vrooj. First, we state a crucial bound that follows immediately form (|2.3H l by l|2.18l l and the fact that ||T"(1 — i?) || < 

Cg(l - gY and i|r|| = 1; namely 

\\Zn{A)\\ < n \\mBiA))\y CAl - ff)l«°°^(^'l (2.37) 

AeA 

Put 

AePo\{N) AeA 

Note that this sum over A £ Pol(A^) can be read as the sum over all A such that exactly one of the A £ A has 
TKookiA) = A^ + 1, and for all other A, THook(^) < N. Moreover, the next largest rHook(^')'^' *= -^ necessarily 
belongs to A. We call this next largest hook time a = THook(^')- ^^ follows that 

Sn<Y m'{B{A))\\# C,(l - g)^^— (-4) Y[{1 + Sa) (2.39) 

AcIn a£A 

where we used that |Hook(A)| = N — max A. Next, observe that Sn is non-decreasing in N, hence one can replace 
Sa on the RHS of l|2.39b by 6n-i- Using g > to control the sum over max(A) and assuming that (1 + 6n-i)£ < 1, 
we get 

6n<^ SUV V (l + <5jv-i)l^l||E^(i?(A))||#<^(l + 5jv-i) (2.40) 

9 roeiN .^, 9 

A C In 
niaxj4 = TO 

where in the last inequality Assumption l2.2l was used. Since 62 < e^ by inspection, we get by induction in N that 
6]\[ < -^—^ for e small enough. Hence eq. | |2.35l l holds. 
To show (|2.36b , we first remark that 

LHSoflESUl < Y II^Jv(^)||Cg(l-g)'"('^)~i (2.41) 

AePo\(N) 



15 



since we get T{1 — R) on the first Tn{A} — 1 factors in the time ordered expression for Z^. Next, we set 

v<{N,n):^ J2 \\Zn{A)\\, v>{N,n) -.^ ^ ll^jv(-4)|| (2.42) 

A G Pol{Ar) A G Pol(Af) 

m{A) < n m{A) > n 

The relations between polymers for different N implies that u> {N, n) — v> [N + k,n + k) for any fc e N. Moreover, 
the summability l|2.35b implies that there is a u < oo such that 

v= lim z;(A^), v{N) -.^ v<{N,n) + v>{N,n) (independently of n) (2.43) 

Given k > 0, we choose N{k) such that \v{N{k)) ~v\<k and hence v<:{N + k,k) < k for any N > N{k) and any 
fc e N. Hence 

Y, \\ZN{A)\\Cg{l~gr^-^^ < vy{N + k,k))Cg{l-gf + v<{N + k,k)) (2.44) 

AePo\{N+k) 

< vCg{l~g)'' + K (2.45) 

As k "> cx), this bound equals k. Since k is arbitrary, this proves (|2.361 . 

D 

2.6.2 Convergence towards the steady state 

From the bounds in | |2.34|| and Lemma [2.51 we obtain 

limsup|lZjv-Zjv'||<limsup ^ \\Zn'{A)R\\+ ^ ||Zjv(^)||Cg(l - g)"(-^) + Cg(l - .g)^ 

^'^°° "^'^"^ ^GPol(7V') ^GPol(JV) 

m{A) <N' - N 

(2.46) 
We will now estimate the two first terms on the RHS of l |2.46|l multiplied by the factor ({N). Note first that for any 
A e Pol(iV') with m{A) < N' ~N, 

N < N' -m{A) + l = Yl ImaxA-minAI + ^ |Hook(^)|, (2.47) 

AgA AeA 

and hence, by property l|1.21l l of the function (,{■), 

({N) < Y[ dist^(yl) X C(|Hook(A)|). (2.48) 

AeA 

Hence 

CiN) Y. \\ZN'iA)\\ < Y ndistc(v4)||i?^(i?(A)|l#C(|Hook(A)|)C,(l-g)l«°°'^(^)l 

A G Pol(Ar') A G Pol(Ar') AeA 

m{A) <N' - N 

< E n c(C,g)l^ldistc(A)||E^(B(A))||#yT^"'°°'^^^' (2.49) 

A G Pol{N') ^6-4 

where we have put (using that C is subexponential, see l ll.23|l ) 



c(C,.9) := sup C{n){^/T^rCg (2.50) 

n>l 

Note that we dropped the restriction that ni{A) < N' — N since it was only necessary for (|2.48t . 
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Consequently, we have derived a bound, (|2.49t , for the first term on the RHS of l|2.46l l (since ||ZAr(^)i?|| < 
||ZAr(y^)||). We will now derive a similar bound for the second term on the RHS of l|2.46ll . 
Instead of (|2.481 , we use here that 

C{N) < CimiA)) X Y[ dist^(A) x C(|Hook(A)|), for A E Pol(iV) (2.51) 

AeA 

and we obtain the same bound as in ll2.491 l, except that this time we get c(C, g)'"^'^^ instead of c(C, g)'"^' because of 
the presence of the term C,{m{A)) in ||2.51| |. 

Next, we show that l|2.49l l (or the analogous bound for the second term on the RHS of (|2.461 l) can be bounded 
by 0{e), for e small enough. To achieve this, we proceed in exactly the same way as in the proof of Lemma [ 
except that here; 



• We include the factor dist|j(A) in the weight ||E'^(i?(T^))||, which is permitted by Assumptions 

• For each set A, there is an additional factor c(C, g)'"^' (or c(C, ff)'^'^^)/ which can be handled by choosing e 
smaller. 



The factor (1 — g) is replaced by ^1 — g, which again forces e to be smaller. 



To conclude the proof of Theorem l2.4i it remains to show that Z,^ '■= linijv/'oo ^7V is of the form I/?™") (1|. This 
follows from the fact that, by the bound (|2.36b , only the terms Zn{A)R contribute to Zoo- CH 

3 Discretization of the physical system 

We explain now how the setup of Section [1] fits into the framework of Section |2l In Section |3?T1 we introduce and 
estimate the Dyson expansion and we present the construction of the operators W{B{A)) from the microscopic 
model. In Sections [3.1.3l and l3.3l we check Assumptions I2.2l and l2.3l starting from Assumptions ll.2l and ll.31 

3.1 Expansions 

To save notation, we present the case where there is only one element in the sum l|1.2|l defining the interaction 
Hamiltonian, i.e. \I\ = 1 and we can write Di = D, ipi — (p. The general case can be treated in essentially the same 
way, we indicate the changes at the end of Section l3.1.2l and in the proof of Lemma r3.2l 
We expand the reduced djrnamics, introduced in Section|2l 

ZnPs = PUiN)U{N^l)...U{l)Pps = TTK[ie-"'''"Hps®Pn'^)e"'''"")] (3.1) 

in a Dyson series. 

3.1.1 Dyson expansion 

We let 

D{t) := e'*^d(i/s)£,e-i*ad(ffs) (3_2) 

and we recall the left and right multiplication operators AfLo(^), MB.i{A) introduced in (|1.12) l. Define the operator 
products 

K{t,k) = Mk,A^D{t2n)) ■ . ■ MkAiDih)), \\K{t,k)\\ < ||i?f" (3.3) 

with t = (ti, . . . , t2n) is an ordered sequence of times < ti < . . . < t2n < t and k = (fci, . . . , k2n) is a sequence in 
{Le, Ri}. Next, recall the correlation functions f^{t) (the labels i, j have been dropped because of the simplification 
\I\ = 1) and define 

r \^f^itr-ts) if ks = Le 

pairings -k (s,r)£7r 



\^f^{tr-ts) if fc^=Ri 
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(recall the convention that s < r in the pairing tt). Finally, we introduce the free S-dynamics 

Wr ■■= e-''^^'^'i(^s) (35) 

Then, we are ready to state the Dyson expansion for Z^; 

Zn = Wn^Y^ I dtK{t,k)G''{t,k) (3.6) 

nSN k o<ti<...<t2„<Niy 

where the term with n = is defined to be 1. This expression can be checked easily by expanding the evolution 
operator e"'*'*'^'^^ ^ in powers of A and writing the expectation values of the field operators in two-point contribu- 
tions by the Wick theorem (this gives rise to the factor G'^(-)). Even though the 'perturbation' H(^-^^ is unbounded, 
it is straightforward to check that the RHS equals the LHS, provided that the sum and integral on the RHS of 03. 6t 
converge absolutely, as will be derived in Section |3.1.2| 

Some explicit intermediary steps of the derivation of | |3.6l l were given in Section [1.7.31 Here we arranged the 
ordered sets of times s_,s^ (cfr. (|1.381 l) into one ordered set of times t and we used the fc-labels to keep track of 
whether a time appeared in the expansion to the left or to the right of the density matrix. 

3.1.2 A formalism for the combinatorics 

The integral over ordered t, together with the sum over left-right specifiers k and pairings, tt, on the set of times, is 
represented as an integral/sum over ordered pairs {{ui, k"), (vi, k^)) with m, Vi E R+ and fc", k" E {Lc, Ri} and 
i = 1, . . . ,n, such that 

Ui <Vi, Ui < . . . < Un (3.7) 

This is done as follows. For any pair {r,s) E vr, we let Ui — tr,Vi — tg and kf ~ kr,k^ = kg where the index 
i ~ 1, . . . ,nis chosen such that the Ui are ordered: -ui < U2 ■ • • < "«■ We represent one pair {{m, kf), {vi^k^)) by the 
symbol Wi and the n-tuple of them by w. We call fi j, with J e IR+, the set of w such that Ui,Vi E J (for arbitrary n), 
and we use the shorthand 

/ dw := ^ ^ du dv x[ui < Vi]x[ui < . . . < u„] (3.8) 

These new coordinates are illustrated in Figure 121 where we relate them to the s, s'-coordinates that were used in 
Section[rz31 

In this new notation, we can write a simple term-by-term bound on the Dyson expansion, using the bound in 
ll331l 

EE / mK{t,k)\\\G\hk)\< j dw\{X^\\D\\f\f\v,~u.)\ (3.9) 

— u<ti <...<r2n<i ^^[o,t] 

The advantage of the last formula is that the sum over pairings is now represented by the integrals over Ui, Vi. One 
can first perform the integrals over Vi, bounding each of them by A^HUlp /^ ds|/'^(s)|, and then the integral over 
Ui gives t"-/n\, i.e. the volume of the n-dimensional simplex. Hence | |3.9|I is bounded by 

e^pUx^D\\hj'\f^{s)\ds] (3.10) 

where the factor 4 comes from the sum over kf, k'f. Since the functions ff^A-) converge to fij{-), uniformly on 
compacts, the Dyson expansion converges absolutely and uniformly in the volume A; 

EE / di\\K{t,k)\\\G^{t,k)\<exp{Xh{\\h\\i+o{\A\°))), asA^Z'' (3.11) 



"6N k o<ti<...<t2„<t 
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Si S2 



■S3 



Si 




Figure 2: Example of a term ("a diagram") in the Dyson expansion. Each bent Hne carries a scalar factor 
/(As) or /(As) where As is the difference between the times at both ends of that line. The dots at the end 
of the line carry operator valued factors D{s), D{s'). In the example, riLc = 4 and nm = 6. We order the 
pairs according to their first time, which is called Ui and the corresponding later time is called u,. The in- 
dex i is fixed by requiring Ui < u^+i. With these (u^, Wi)-coordinates, the above diagram would correspond to 
{ui,vi) = (si,S2),(u2,W2) = (s'2,s4),(u3,u3) = (si , S3), (u4, W4) = (sg, Sg), (^5, U5) = (S4, S5) The left/right labels 
kf, fc" would then be ki,k2,k^,k^ = Lc and fc", k2,k^,kl^,kl,k^ = Ri. That is, whenever a time is on the upper 
line, it has the label Le and, whenever it is on the lower line, it has the label Ri. 



where the positive function h{t) was defined in (|1.151 for a general index set I. It follows that (here we indicate 
explicitly the A-dependence of Z^) 



lim^ Z^=^WnJ2 I '^^-Yl ^(^' ^)'^(^' ^) 



(3.12) 



0<ti<...<t2„<i^A^ 



where G{t,k) = liui^y^^d G^{t,k) is obtained by simply replacing /^(•) by /(•). We have now proven Lemma im 
up to the simplification that \I\ = 1. However, one can trivially check that in the case \I\ > 1, the bound l|3.11t 
remains true. Finally, we introduce the correlation function 



Giw) 



n 



\^J{v-u) if fc" = Le 
\^f{v-u) if fc" = Ri 



(3.13) 



where the product runs over the pairs ((w, fc"), (w, k"")) in w_. 

3.1.3 Connected correlations and the Dyson series 

To relate the previous sections to the setup in Section |2l we need to discretize time and express the operators 
W{B{A)) in terms of the Dyson series. 

Recall that /at — {1, ... , N} is the set of macroscopic times. To a set A <z In oi macroscopic times, we associate 
the domain of microscopic times 

Dom,(A) = IJ [v{t - 1),vt] (3.14) 

reA 

A set of pairs w_ G ^Y)oai^(A) determines a graph Qa{w) on A by the following prescription: the vertices t < t' are 
connected by an edge if and only if there is pair w = ((u, k"), {v, k")) in w such that 



u e Domy(r) = [i^(r — 1), i^t] and v E Dom^(r') = [i^(t' — 1), i^t'] 



(3.15) 



We write Supp(tjy^ (w) ) for the set of non-isolated vertices of Ga {ul), i-e. the vertices that have at least one connection 
to another vertex, (this set Supp{Qa{w)) is obviously a subset of A). 
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Now, our task is to connect the Dyson expansion in sets of pairs with the discrete-time expansion that ^vas 
introduced in Section|2l The connection is via 



(g) W-r E(B(A)) (^ Wr-l 



(3.16) 



dw^Ir(t.„) [Mfe,„(iD(i2„))] . . .!.(*,) [Mfc,(ii^(ii))] G{w) 



(3.17) 



Supp(eA(2«)) = A 



where Ir is the embedding that maps operators on S^ into operators on S^r, and T{t) assigns to each microscopic 
time t the right "interval" t (that is t = r(t) o t e Domy(T)). The coordinates ii, fci on the RHS are determined 
imphcitly by the w_ as explained in Section [3.1.21 The product of correlation functions G(w) has been defined in 

The above formula follows immediately from the definition of i? ( • ) . Indeed, sets of pairs u; such that Supp ((J^i ( w) ) 
is strictly smaller than A, say r e ^\Supp(tJ^(w)), can be resummed to giver(r), and hence they do not contribute 
to B(t). Indeed; 



(3.18) 



W^rT{T)W,^X = I dw I,(t,„) [Mk„AiD{t2n))] • • -^rit,) [Mfe, (ii?(tl))] G{w) 

where now all T{tj) = t and hence all embeddings I^(^t j are into Mr- 
Next, we state a useful formula for the "truncated correlation functions". 

Lemma 3.1. Assume the simplification \I\ — 1, then 



Qg W-r E-{B{A)) Q^ Wr-l 
VreA / \tGA / 



j dwlrit,„) [Mk,JiD{t2n))] ■ . -Irit^) [Mfe, (iD(ii ))] G{w) 



(3.19) 



W e flDom,AA) 

Qa (w) connected 

See Figure|3]for two examples of w that contribute to W^{B{A)) for a given A. 

Proof. Let 

F{w):=K{t,k)G{w) (3.20) 

where, again, the t, k are determined implicitly by w, as explained in Section l3.1.2l Obviously, we have the decom- 
position. 



/ 



dwF{w) 



(3.21) 



m S f^Dom,AA') 
Supp(e^,(«,)) = A' 



E 

A partitions of A' 



n 

AeA 



dWA 



Qa {w_a ) connected 



F{ U w^) 



(3.22) 



where U w_j^ is well-defined since w.a ^^^ sets (of pairs of times). Since the connected correlation functions are 
uniquely defined by (|2.141 , the claim of the lemma follows. D 
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I- 



A^ 




Figure 3: We illustrate the discretization of the perturbation series. The numbers 0, 1, . . . , A^ refer to macro- 
scopic times measured in units of v. The top drawing shows (by the graph under the horizontal axis) the set 
A — {1, 6, 7, 9}. This set in fact refers to the time-intervals that end in these points, i.e. the intervals that are marked 
by horizontal lines. In the two bottom pictures, we have drawn two examples of diagrams (without indicating the 
fc-labels) that induce a connected graph on A and hence contribute to W{B{A)). 



3.2 Estimates 

In this Section, we verify Assumption 12.21 of Section |2l The following estimate is an immediate consequence of 
Lemma |3A] and the bounds at the end of Section |3.1.1| 

Lemma 3.2. 

||E"(S(v4))||# <e^'ll'*"il°°'"''('^)l / Aw i\[\^h{vi~ui)\ x[w spans Aminimally] (3.23) 

where the statement "w spans A minimally" means that Qa{w) is connected and that no pair can be dropped from w without 
losing this property. In particular, this implies that Qa{w) is a spanning tree on A. 



Proof. Let F{-] be as in l|3l20t , then 

/ Aw\\F{w)\\ < 



I 



m S f2Dom„(A) 

Qa {w) connected 



M' e f^Dom„{A) 

u' spans A minimally 



du;' / dw"||i^(«;'U2z;")|| 

m"GODom„(/l) 



(3.24) 



This appealing estimate was the main motivation for encoding the pairings tt in the elements w_. 

To realize why it holds true, choose a spanning tree ,'7 for the connected graph Qa {w) and then pick a minimal 
subset w' of the pairs in w_ such that Ga{w!) = ^ ■ Since, in general, this can be done in a nonunique way, the 
integrals on the RHS contain the same w_ more than once, and the inequality is strict unless F is concentrated on 
minimally spanning w_. Starting from this inequality, we now use the bound l|3.11t to perform the integral over w", 
using also that ||-F(u))|| < nr=i ^^H^i — "O- ^ 

Remark. This lemma holds in the general case \I\ > 1. To see this one can first extend the setup slightly, by making 
wa set of pairs (u, k'^,a"';v,k'",a") where a", a" range over the elements inX. Then, the reasoning goes through in 
exactly the same way and at the end the sum over a" , a" can be performed (note indeed that the function h includes 
a double sum over I, see ||1.15|| ) such that the bound Il3.23b remains true with the original definition of w. 
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We now prove that Lemma 13.21 implies Assumption 12.21 Recall the decay time t^ associated to the Markov 
approximation by Assumption [T31 

Proposition 3.3. Set the scale factor v = \^^£ where A is the coupling strength and £ £ [i^, 21^]. Fix also a nondecreasing 
function C(-) satisfying (|1.21|l and such that J dt(^{t)h{t) < oo. Then, for A small enough there exists e — e(A) with e(A) — > 
as |A| — >■ such that 

sup J2 e-l^ldistc(A)||E^(B(A))||# < 1 (3.25) 

^°^'" AcIm 

max(A) = TQ 

uniformly for i e [i^, 21^] and N. 

Proof. Each w that spans A minimally determines a spanning tree on A. Hence we can reorganise the bound (|3.231 
by first integrating all w that determine the same tree. This amounts to integrate, for each edge of the tree, all pairs 
{u, v) that determine this edge. Hence we arrive at the bound 

span. trees S^ on A {t,t'}^£{^) 

where S{^) is the set of edges of the tree J^, 

e{T,T') := /" dv f dvX'^h{v-u), if r < r' (3.27) 

Ju[t,t + 1\ Ji/[r',T'+i\ 

and e(r', r) :— e(T, t') for r' < r. To bound e(r, r'), we use 

i/(r'-r+l) 



Kt^t') < X^v / ds/i(s), for t' - t > 1 



i/(r''— r — 1) 
2y 



e(r,r+l) < AM dssft,(s) 
Jo 



(3.28) 



(3.29) 



Using the integrability assumption / dtC^(t)h{t) < oo, the inequality ({\t — t'\) < ({\t — t'\ — l)C(l) from l|1.21|l , and 
the fact that ( is nondecreasing, we easily derive that 

Y^ C(k-^'l)e(r,r') <ei, with ei = ei(A) ^ as |A| ^ (3.30) 

r'eN\{T} 

uniformly ini E [t^, 2ii]. The key observation here is that for any integrable function /(•), 

1 /•* 



, dssfis) — ^ (3.31) 

t Jq t->00 

which follows by the dominated convergence theorem applied to the sequence of functions Ft{s) := x[s < i]f /(s), 
since V.s : limf/ioo Ft{s) = 0. Note that our bound I I3.30D is not very ambitious in that we did not try to use the fact 
that C, increases already on the microscopic scale. 

We now turn to the sum in | |3.26l l. Recall the distance functiondist^; ( A) introduced in Section [231 Since C is an 
increasing function, the inequality 

distc(A)< n C(|r'-r|) (3.32) 

{r,r'}G£(.T) 

holds for any spanning tree ,3^ on A. Here is a procedure for checking | |3.32| |: Let En^ be the set of edges of the 
linear tree, i.e. £iin = {{tz, ^i+i}, « = 1, • ■ • , |^| — 1} where ti, . . . , t|^| are the time-ordered elements of A. Then we 
need to prove 

n '^e < n ^-^ = di£i.r)), d{,,,,} = C(|t' - r|) (3.33) 
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Choose a leaf r of the tree S^ and let r' be the unique vertex that shares an edge (that we call e^) with r. If e^- e E\\^, 
then (|3.33b is equivalent to the claim where we replace fun by E\\^ \ e^ and S' by the tree with the vertex r and the 
edge er removed. If Cr ^ ^\m then consider the tree S' where the edge e^ is removed and the edge er ~ {t, t} is 
added, where f is chosen such that f — r has the same sign as r' — t and Bt G f lin- Clearly, (i(£(^)) < d{£{,!7)) and 
hence it suffices to prove (|3.331 with £^ replaced by i^. We can continue by induction. 

Using that A^ |Dom^(A) | < 2tL | A| in the exponent on the RHS of (|T26t , the LHS of (|3l25l is bounded by 

sup ^ ,-|A|eCo„.t|A| ^ ]-j ^(|^,_^|)g(^^^,) (334) 

'^"^'^ A C N : A 3 TO, |A| > 2 span.trees5^onA {r,T'}e£(.S^) 

(Note that we also replaced the constraint max A = tq hy A 3 tq). The sum over A 9 tq and spanning trees on A is 
of course equivalent to a sum over trees containing tq. The control of such sums is at the core of cluster expansions, 
see e.g. [29J. We state a simple result. Lemma [3.41 proven in a more general setting in |29| . D 

Lemma 3.4. Let v{-,-) he a symmetric and positive function on N x N such that 

sup y^ v{t,t') <-e~'^, for somen >0 (3.35) 

"^ t'£N\t 

Then 

sup ^ J]^ v{t,t')<k (3.36) 

^" trees^onN (T,r')e£(5^) 

1< |5^| < oo,5^9to 

We apply this lemma to our case with w(r, r') = e^^Cd''"' — '''|)e(T, r') and k = gconst^ f^^ ^ sufficiently small. 
It follows that l|3.34l l is bounded hy 0{^). This means that one can take e — > as A — > 0, since ei (A) -^ 0. 

3.3 The dissipativity condition on T 



In this section, we check Assumption l2.3l Observe first that the reduced dynamics T = Tx^i depends on both A and 
£ (via the scale factor i/). We quote a celebrated result (originally in ||6l under slightly more stringent assumptions, 
see irniTOl for the result as quoted here) 



Theorem 3.5 (weak coupling limit). Assume that Assumption \1.2\ holds. Then, for any (.max < oo, we have 



lim sup 
where L is the generator given in Section^ 



gi£A-^ad(Hs)y^^_gfL 



(3.37) 



Remark In fact, the proof of Theorem l3.5l is rather straightforward if one uses the framework of the present paper. 
In l|3.18ll , one shows that the contribution of those w_ that contain two pairs (wi, u^), (wj, Wj), J < j such that Uj < Vi 
(i.e. the pairs are 'entangled') vanishes in the limit A ^ 0. This follows from the same simple calculation as done 
in the proof of Proposition 13 .31 to show that ^^, e(r, r') vanishes as A — > 0. Then, all what remains is to prove that, 
the 'ladder diagrams' Wj i.e. those for which Vi < u^+i, sum up to produce the semigroup c^ *^, in the interaction 
picture and up to an error that vanishes as t = A^^t, A — > 0. 

By Assumption II. 3[ we know that the semigroup c'^ is exponentially ergodic (see l|1.17|l '). Since the free dy- 
namics Wi := e^'^^ ad(ifs) jg gj-( isometry on ^i (J^s), we deduce from (|1.17t that 

WW.e''^ - Wi\pl){l\ II < e-'/'^ fori > iz, (3.38) 
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Since L commutes with ad(iJs)/ we also have Wep^ = pg. Hence, to check Assumption 12.31 we merely need to 
choose A small enough such that 

Tx,i - Wee'"^ (3.39) 

is a sufficiently small perturbation of W^e^^ . In that case, spectral perturbation theory of isolated eigenvalues 
applies and we find that 7\ £ has a unique maximal eigenvalue (which, by unitarity, is necessarily equal to 1) 
corresponding to the eigenvector |ps (^))(1| where pg — Pg (^) = o(|A|°), as A ^ 0. The appearance of (1| and the 
fact that Pg {€) is a density matrix, are consequences of the fact that the dynamics preserves positivity and the trace. 
Moreover, 

ll(TA,.r-ipi'w)(iill<(i-.g(^)r (3.40) 

where 

1 - .g(£) = e-^A- + o(|A|") = exp {--— ^— ^}, A ^ (3.41) 

tL+0(|A|") 

In the above estimates, the error term o(| A|") depends on I as well; in general, a smaller A is needed when I grows. 
We restrict ^ e [ti, 2ti], which allows us to make all estimates uniform in i. 

The conclusion is that, for A small enough. Assumption 12.31 holds for T = T\^g, with £ e [1^,21^], g(£) as in 
^M and Cg = 1. 



4 Proof of Theorem ILH and Proposition [L5 



Finally, we weld the abstract results of Section |2] to the setup presented in Section [TJ using the estimates of Section 
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4.1 Theorem [M 



First, we choose A small enough such that 

• Assumption 12.31 holds with T = T\j for I E [II, 2ii] and such that g can be chosen uniformly in £. The 
possibility of doing this was sketched above in Section l33l 

• Proposition 13.31 applies with e(A) small enough such that Assumption 12.21 holds and Theorem l2.4l applies. 
Let us abbreviate the reduced dynamics Vt := Pe~'*^'^^^^'P, then Theorem |23] yields that 

Vmx-^- --^ IprWXll- (4-1) 



N- 



Moreover, we deduce that p™^ {£) is continuous in £. This follows from the absolute summability of the series for 
Zco (by l l2.35|l ), uniformly in £, and the continuity of W{B{A)) and T as functions of i (which follows easily from 
bounds on the Dyson expansion). 

Next, we argue that p™'"{£) does not depend on £. Indeed for £i, £2 such that £1/^2 is rational, we can pick a 
subsequence of (Vjv^iA-^)^ that is identical to a subsequence of {Vni^x-2)n. Hence the limits of both sequences 
coincide. Since £ —>■ pg"" {£) is continuous, it is necessarily constant. This ends the proof of Theorem ll.4l D 



4.2 Proposition IT3 



We start from (|2.28| l where we put Cg = 1 (as argued in Section l331 this is possible) 

\\V^,x-2-\pr){l\\\<{l-9{e)f + 0{e)^y e^O (4.2) 

Since any t > X~^iL can be written as t = X^^N£ for some A^ > and £ e [t^, 2tL], we obtain 

11^* - IPDdlll < -p{-i7^^} + ^W^))c(a4m' ' - « (4.3) 

where we used ( |3.41t . After using that e ^ as A ^^ 0, and C,{\^t/£) > ({X^t/{2iL)), we obtain Proposition [LS] D 
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